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Abstract — We consider the wireless Rayleigh fading relay 
channel with finite complex input constellations. Assuming global 
knowledge of the channel state information and perfect syn- 
chronization, upper and lower bounds on the achievable rate, 
for the full-duplex relay, as well as the more practical half- 
duplex relay (in which the relay cannot transmit and receive 
simultaneously), are studied. Assuming the power constraint at 
the source node and the relay node to be equal, the gain in 
rate offered by the use of relay over the direct transmission 
(without the relay) is investigated. It is shown that for the case 
of finite complex input constellations, the relay gain attains the 
maximum at a particular SNR and at higher SNRs the relay 
gain tends to become zero. Since practical schemes always use 
finite complex input constellation, the above result means that 
the relay offers maximum advantage over the direct transmission 
when we operate at a particular SNR and offers no advantage 
at very high SNRs. This is contrary to the results already known 
for the relay channel with Gaussian input alphabet. 

I. INTRODUCTION 

The study of the wireless relay channel has attracted a lot 
of attention in recent times. Achievable rates and the upper 
bounds on the capacity for the full duplex (FD) discrete mem- 
oryless relay channel (DMC) were studied in [ 1 1. Motivated by 
the practical constraints involved in the construction of relays 
which can transmit and receive simultaneously, the bounds 
on the capacity for the so-called cheap relay networks were 
derived 0. We will refer to such a relay channel as the half- 
duplex (HD) relay channel in this paper. These results were 
extended to the Gaussian relay channel with continuous input 
and output alphabets Q. 

For a single-user Gaussian channel, if the input is con- 
strained to come from a finite complex constellation with the 
points used with uniform distribution, the mutual information 
between the input and the output of the channel is called 
the Constellation Constrained Capacity (CCC) in |4). Similar 
notion of CCC was studied for the Gaussian-MAC channel in 
1 5 1 and for the Gaussian Interference Channel in [6 1 with finite 
complex input constellations. In this paper we study bounds 
(both upper and lower bounds) on the achievable rate of the 
relay channel with finite complex input constellations. 

The lower bound on the achievable rate of the relay channel 
is defined as be the achievable rate of the decode and forward 
scheme Q), when the input symbols take values from a finite 
constellation. Similarly, the cut-set bound given by the max- 
flow min cut theorem Q, when the input symbols take values 



from a finite constellation, is referred to as the upper bound 
on the achievable rate. If the expressions for the bounds 
involve maximization with the respect to the joint distribution 
of random variables, it is assumed that the random variables 
are independent and uniformly distributed. For the fading 
relay channel, since the mutual information expressions which 
appear in the bounds are functions of the fading coefficients, 
expectation with respect to the fading coefficients is taken. 

Results from 0, where the capacity bounds and optimal 
power allocation strategies for the wireless FD and HD relay 
channel were studied, reveal that for the Gaussian input 
alphabet, the gain in rate offered by the relay channel over 
the direct transmission without the relay node, referred to as 
the relay gain, is significant and in particular at high SNR, 
this gain becomes a constant. 

The main contributions of the paper are as follows. 

« We present lower and upper bounds on the achievable rate 
of the relay channel for both the FD and HD scenarios. 

• It is shown that for the case of finite input constellations, 
the relay gain attains the maximum at a particular SNR 
and at higher SNRs the relay gain tends to become zero. 
Since practical schemes always use finite input constel- 
lation, this result means that the relay offers maximum 
advantage over the direct transmission when we operate 
at a particular SNR. 

• The variation of the relay gain with the duty cycle of the 
HD relay channel is studied. 

The organization of the paper is as follows. Section II 
describes the FD and the HD relay channel models. In Section 
III, we briefly discuss the already known capacity bounds for 
the relay channel and then we derive these bounds for input 
constellation constrained FD and HD fading relay channels. 
In Section IV, the bounds and the relay gain with 4-QAM 
as the input constellation are computed. A comparison of 
these results with already known results for the Gaussian input 
alphabet is made. Also, the variation of the relay gain with the 
duty cycle of the HD relay channel is studied. 

Notations: For a set S, \S\ denotes the cardinality of S. 
For a random variable X s which takes value from the set S, 
we use x St i to represent the i-th element of S. E Z [Y] denotes 
the expectation of Y with respect to the random variable z. 
Throughout, log refers to and C(a) denotes log(l + a). 




Fig. 1. The Relay Channel 

II. CHANNEL MODEL 

We consider the relay channel shown in Fig. Q] consisting 
of the source node S, the relay node R and the destination 
node D. We focus on the Rayleigh fading scenario and let 
h rs = c rs e^ r " , hd s = Cd s e^ da and h dr - c dr e^ dr denote 
the zero mean complex Gaussian fading coefficients associated 
with the source-relay, source-destination and relay-destination 
links respectively, with the corresponding variances denoted by 
a rs , ads and Odr- It is assumed that the power constraint at the 
source and the relay to be equal, denoted by P. Throughout, it 
is assumed that the global information about the channel state 
information (CSI) is available, i.e., all three nodes know the 
instantaneous values of h rs , h ds , h dr - It is also assumed that 
all the nodes have perfect timing and carrier synchronization. 




Fig. 4. HD Relay Channel - Phase 2 

The following two scenarios are considered: full-duplex 
(FD) relay channel (Fig. |2), in which the relay node can 
receive and transmit simultaneously and the half-duplex (HD) 



relay channel, in which the relay receives during a fraction of 
the time aT, and transmits for the duration (1 — a)T, where T 
denotes the frame duration. In this paper we refer to a as the 
duty cycle of the relay. The receive and the transmit phases of 
the relay are referred to as Phase 1 (Fig. [3} and Phase 2 (Fig. 
|4]i respectively. Since we have assumed that CSI is available 
at all the three nodes, the phases can be cancelled out at the 
appropriate places. The phase 4>ds due to the S-D link can be 
cancelled out at S and the resulting phase <fi rs ~ <fids at R due 
to the S-R link can be cancelled out at R. The phase <fidr due 
to the R-D link can be cancelled out at R. Hence we can take 

4>rs=4>ds=4>dr=Q- 

For the FD relay the received signal at R and D are given 

by, 

Y r c rs X s -\- z r , 

Yd = Cd s X s + Cd r X r + Zd- 

where X s denotes the symbol transmitted by S, z r and z s are 
independent complex Gaussian random variables with mean 
and variance 1/2 per dimension, denoted by CG(0, 1). 

For the HD relay, the received signal at R during Phase 1 
is given by, 

The received signal at D during Phase 1 and phase 2 are given 

by, 

=c ds Xl+z d ; 

Yj = C ds X 2 s + C dr X r + Z 2 d . 

where X} and X 2 denote the symbols transmitted by S during 
Phase 1 and Phase 2 respectively, X r denotes the symbol 
transmitted by R during Phase 2. z\, z\ and z r are independent 
and CG(0, 1). 

III. BOUNDS ON THE ACHIEVABLE RATE 

The bounds on the achievable rate for fading relay channel 
with BPSK input were plotted in Q OH. In this section, the 
bounds for the fading relay channel with an arbitrary finite 
complex input constellation are derived. To derive the bounds, 
we follow a similar procedure as in |]5] |0. 

A. FD relay 

The lower and upper bounds on the capacity for the discrete 
memoryless FD relay channel JT] Q are given by, 

Clfd =maxp X3:Xr min{I(X s ,X r ;Y d ),I(X s ;Y r \X r )}; (1) 
Cufd =max FXs:Xr min{I(X s ,X r ;Y d )],I(X s ;Yr,Y d \X r )}. (2) 

Throughout, it is assumed that X s and X r are independent and 
hence the max operation in (03 and (f2|l is left out in the rest of 
the paper. For the FD relay with Gaussian alphabet, assuming 
that X s and X r are independent, the lower and upper bounds 
on the capacity are given by JS) 



Cf FD = min {E Cds , Cdr [C ((c 2 da + c 2 dr )P)] , E Cr3 [C («£.P)] } ; 
Cg FD = ram{E Cj „ eir [C ((c ds + c 2 dr )P)], 

Ec rs ,c ds [C((c 2 ra + 4 s )P)]}. 



To compute the bounds on the achievable rate, it is assumed 
that X s and X r are independent and uniformly distributed, 
with the same power P, X s G S s , X r G S r , with \S S \ = M s 
and \S r \ = M r . 

The achievable rate bounds of the FD relay channel with 
finite input constellation are given by, 



Clfd = min{Ec dg ,c dr [Ri],E CrB [R2]}; 



Cufd = min{Ec d3 ,c dr .[Ri],E CrB , Cds [R3]}. 



(3) 
(4) 



Expressions for i?i,i?2 and R3 are given in $1% , (l22l 
and d25l l respectively along with their explicit calculations in 
Appendix lAl The plot of the achievable rate bounds Vs P, for 
<Tds = -10 dB, a rs = 2 dB and a^r = 12 dB, for the FD relay 
is shown in Fig. [5] (in the next page). In Fig. [8] a similar plot 
is shown for <7d s =-10 dB, <j rs =l2 dB and Udr=2 dB. 

B. HD relay 

For the discrete memoryless HD relay channel with duty 
cycle a, the upper and lower bounds on the capacity are given 

by da, 



C LH d = maxp , 



x^xi ,x r 



min{al (X a ;Y r 



c, 



uhd = maxp 



Xi ,xi ,x 



+ (l-a)I(X 2 ;Y d 2 \X r ), 
a I(Xl- Yj) + (1 - a)I(X 2 , X r ; y|)}; (5) 
min{aI(Xl;Y r , Y d ) 

+ (l-a)I(X 2 ;Yi\X r ), 
al (Xl ; Y} ) ] + ( 1 - a) I (X 2 , X r ; Y d 2 ) } . (6) 

Similar to the FD relay, it is assumed that X^, X^ and X r 
are independent and hence the max operation in (0 and (0 
are left out. 

For the HD relay with Gaussian alphabet, assuming that 
X] , X"l and X r are independent, the lower and upper bounds 
on the capacity are given by O 0, 



CShd = min {aE Crs [C (c 2 rs P)] + (1 - a)E Cds [C (c 2 da P) 

r 2 P 

c dr jr 



aE Cda [C (c 2 da P)] + (l-a)E Cd3 , Cdr 



c Cd _ s p + 



Cuhd = min {aE Cds , Crs [C ((c 2 ds + c 2 , s )P)} 

+ (1 - a)E Cds [C (cLP)] 
r 2 P ' 



(1 - a)E c 



C[c' da P + 



1 



+aE Cds [C(c 2 ds P)]} 



Note that we have assumed the same power constraint P for 
both S and R. S uses the same power P in both Phase 1 and 
Phase 2. Since the relay transmits only for a fraction (1-a) 



of the total time, the transmit power P of the relay has to be 
scaled by a factor l/(l-a). 

To compute the achievable rate bounds with finite input 
constellation, assume X] G Si, with \Sl\ = M Sl , X^ G S^, 
with \S% I = M S2 and X r G S r , with \S r \ = M r . 

Define R A = I(Xl;Y r ), R 5 = I{X 2 s -Y^\X r ), R 6 = 
I(Xl;Yi), R 7 = I(XlX r ;Yl) and i? 8 = I(Xl;Y r , Y^). 

The achievable rate bounds of the HD relay channel with 
finite input constellation are given by, 

Clhd = min{aE Cr3 [R 4 ] + (1 - a)E Cds [R 5 ], 

a E Cds [7? 6 ] + (l-a)E Cd „ Cdr [ J R 7 ]}; (7) 

Cuhd = min{aE Crs , Cds [R 8 ] + (1 - a)E Cds [R 5 ], 

«E Cds [i? 6 ] + (l-a)E Cds , Cdi .[ J R 7 ]}. (8) 

Expressions for R4 , R 5 , R e , R? and R s are given in (|35l l- 
( |39l > along with their explicit calculations in Appendix iBl The 
plot of the achievable rate bounds Vs P for the HD relay, for 
Cds=-10 dB, <r rs =2 dB, (jrf. r =12 dB and a=0.5 is shown in 
Fig- 121 (in the next page). In Fig. [8] a similar plot is shown for 
fds=-10 dB, <r rs =l2 dB, (Jdr=2 dB and a=0.5. 

C. Direct transmission 

Our goal is to compare the bounds for FD and HD relay 
with the capacity of direct transmission without the relay and 
in order to make the comparison fair, we assume that the power 
constraint at S is 2P. For direct transmission with Gaussian 
input alphabet, the capacity is given by, 



^DIR 



= E Cds [C(2c 2 ds P)} 



and to compute the CCC with finite input constellation, 
assume X s G S s , \S S \ = M s and define Rg = I(X s ;Yd). Rg 
can be computed and is given by (0 (in the next page). The 
CCC for the direct transmission without the relay, with finite 
input constellation is given by, 

Cgf R = E[R 9 ]. (10) 

Fig. [5] and Fig. [7] show a comparison of C^,f R with the 
achievable rate bounds of the FD and the HD relay channels. 

IV. RELAY GAIN COMPUTATION 

Definition 1: Relay gain is defined as the difference be- 
tween the lower bound on the achievable rate in the presence 
of the relay and the capacity of the direct transmission without 
the relay. 

In this section we compute the achievable rate bounds and 
the relay gain for the FD and HD relay channels. The capacity 
bounds and relay gain for the fading FD and HD relay channels 
with Gaussian input alphabet were studied in [8 1. It is assumed 
that the finite input constellation used at S and R are the same, 
i.e., M s = M r = M for the FD relay and M Sl = M S2 = M r 
= M for the HD relay. 

The expressions for the capacity bounds, given in Appendix 
lAl and Appendix [B] can be evaluated using Monte-Carlo 
simulations. Fig. [5] and Fig. [6] show the Rate Vs P plot for 



the FD relay channel, showing the lower and upper bounds 
on the capacity with Gaussian and 4-QAM input alphabets, as 
well as the capacity of the direct transmission with Gaussian 
and 4-QAM input alphabets, for two different sets of cr<j s , &rs 
and <Tdr- In Fig. |7]and[8] similar plots are shown for the HD 
relay with a=0.5. 

From Fig. [5] Fig. [6] Fig. [7] and Fig. [8] under the given 
assumptions, we observe that the upper and lower bounds on 
the capacity are very close, not only for the Gaussian input 
alphabet, consistent with the results in (8), but also for the 
finite input constellation (4-QAM). At low values of P, the 
rate achievable with the use of finite input constellation is 
nearly the same as the rate achievable using Gaussian input 
alphabet. 




Gaussian Alphabet 



—Relay Lower Bound-4 QAM 
■•■Relay Upper BounrJ-4 QAM 
Capacity Direct Tx.-4 QAM 
Relay Lower Bound-Gaussian Alphabet _ 
■▼■Relay Upper Bound-Gaussian Alphabet 
■■-Capacity Direct Tx -Gaussian Alphabet 



10 15 
PinrJB 



25 



- ' ' /I ' /I 

// / Gaussian Alphabet - 
/'* /" 4 QAM 

/' c - - 


• * .--i^ — * 




—Relay Lower Bound-4 QAM 
■•-Relay Upper Bound-4 QAM 
Capacity Direct Tx.-4 QAM 
""- Relay Lower Bound-Gaussian Alphabet 
-T- Relay Upper Bound-Gaussian Alphabet 
■■-Capacity Direct Tx.-Gaussian Alphabet 



Fig. 7. Rate Vs P for HD relay with <r ds =-10 dB, cr rs =2 dB, a dr =i2 dB, 
a=0.5 
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Fig. 5. Rate Vs P for FD relay with a dB =-W dB, oy s =2 dB, a dr =12 dB 
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Fig. 6. Rate Vs P for FD relay with cr ds =-10 dB, oy 3 =12 dB, <r dr =2 dB 



Fig. 8. Rate Vs P for HD relay with o- ds =-10 dB, o> s =12 dB, a dr =2 dB, 
a=0.5 



A. Variation of relay gain with P 

The relay gain, which is the difference between the curves 
corresponding to the achievable rate in presence of relay and 
the capacity of the direct transmission, in Fig. [5] Fig. [6] Fig. [7] 
and Fig. [8] is plotted in Fig.|9] Fig.[lO] Fig. [U] and Fig. [12] As 
observed in JSJ, with Gaussian input alphabet, the relay gain 
increases with P initially and at high values of P, becomes 
a constant. Expressions for the asymptotic relay gain (relay 
gain as P ^ oo) for the FD and HD relay with Gaussian 
input alphabet were derived in (8). 

Contrarily, with finite input constellation the asymptotic 
relay gain is always zero. Assuming that the finite constellation 
used at S and R are of the same cardinality A/, it can be easily 



seen from the equation sets {©.dHJ,©}, {@IM38),l|7}} 
and {d9li.([T0li} that as P — > oo , the achievable rate with the FD 
and HD relays as well as the capacity of the direct transmission 
become equal to log(M) and hence the asymptotic relay gain 
is zero. As seen in Fig. [9] Fig. [10] Fig. Q~TJ and Fig. [12] the 
relay gain with finite input constellation, attains a peak value 
at a particular value of P and becomes zero at high values of 
P. The same argument holds good for arbitrary values of ads, 
a rs and adr- This peak value of the relay gain is referred to 
as the maximum relay gain. 

It is interesting to note that similar conclusions were ob- 
served for the secrecy capacity of the wire-tap channel with 
Gaussian and M-PAM input alphabets in (Til and with QAM, 
PSK input alphabets in fl2ll . 



Fig. 11. Relay Gain Vs P for HD relay with a ds 
dB, a=0.5 
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Fig. 9. Relay Gain Vs P for FD relay with a ds =-W dB, a ra =2 dB, a dr =12 Fig 12 . Re i ay Gain Vs P for HD relay with o" ds =-10 dB, oy s =12 dB, cr dr =2 
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Fig. 10. Relay Gain Vs P for FD relay with cr da =-K) dB, oy s =12 dB, o dr =2 
dB 



B. Variation of Relay gain for HD relay with a 

Fig. [13] shows the Relay gain Vs P plots for four different 
values of a, for 4-QAM input alphabet, with <Jds = -10 dB, 
a rs = 2 dB and a^r = 12 dB. As seen from Fig. [13] the 
maximum relay gain increases when a increases from 0.3 to 
0.7. When a increases from 0.7 to 0.75, the maximum relay 
gain decreases. This means that the maximum relay gain as 
a function of a, attains a peak for some a between 0.7 and 



0.75. Similarly, from Fig. [14] we can see that for ads = -10 
dB, a rs = 12 dB and Odr = 2 dB, the maximum relay gain as 
a function of a, attains a peak for some a between 0.3 and 
0.5. 

Fig. Q3] shows the Relay Gain Vs P plots for four different 
values of a, for Gaussian input alphabet, with ad s = -10 dB, 
a rs = 2 dB and a^r = 12 dB. As seen from Fig. Q3] the 
asymptotic relay gain increases when a increases from 0.3 to 



■►■Duty Cycle = 


0,3 


■►Duty Cycle = 


0.5 


■♦■Duty Cycle = 


0.7 


■ Duty Cycle = 


0.75 



10 15 
P in dB 



Fig. 13. Relay Gain for HD relay for 4-QAM with a ds =-iO dB, oy s =2 dB, 
a dr =\2 dB for four different values of a 
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Fig. 14. Relay Gain for HD relay for 4-QAM with cr ds =-10 dB, a rs =12 
dB, u dr =2 dB for four different values of a 
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Fig. 15. Relay Gain for HD relay for Gaussian alphabet with ct £ j s =-10 dB, 
<T rs =2 dB, <r t j r =12 dB for four different values of a 



0.7. When a increases from 0.7 to 0.75, the asymptotic relay 
gain decreases. The asymptotic relay gain as a function of a, 
attains a peak for some a between 0.7 and 0.75. Similarly, 
from Fig. [16] we can see that for ads = -10 dB, a rs = 12 dB 
and <Tdr = 2 dB, the asymptotic relay gain as a function of a, 
attains a peak for some a between 0.3 and 0.5. 

From Fig. [13] Fig. [14] Fig. [15] and Fig. [16] it is very clear 
that the variation of the relay gain with P and a is not the 




same with Gaussian alphabet and 4-QAM. The fact that with 
4-QAM, the maximum advantage over the direct transmission 
is obtained at a particular value of P and a cannot be seen 
from the relay gain plots for Gaussian input alphabet (Fig. [T31 
and Fig. [16}. 

V. DISCUSSION 

It is shown that the variation of relay gain as a function 
of SNR with finite input constellation is different from the 
variation with continuous Gaussian input assumption. It will 
be interesting to study how the Relay gain Vs SNR curves vary 
with increasing size of the finite input constellation. Further it 
will be interesting to study the channel conditions under which 
the relay gain with finite input constellation is negative, i.e., the 
use of relay offers no advantage over the direct transmission 
without the relay. 
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Appendix A 

Calculation of bounds for the FD relay channel 

WITH FINITE INPUT CONSTELLATIONS 

Define R x = I(X 8 , X r ;Y d ), R 2 = I(X s ;Y r \X r ), and 
R 3 = I{X s -Y ri Y d \X r ). 

1) Calculation of R\: For the FD relay, I(X S , X r ;Yd) is 
given by, 

I(X s ,X r ;Y d ) = h(Y d ) - h(Y d \X s ,X r ), 

where h(Yi) can be calculated using the density P(Y d ) given 
in (fl2l > and is given by (TlTb and 

M,-ijf r -i 

h(Y d \x s ,x r ) = E h{ y*\ x ° = *s,i,x T = x r ,j), 

where h(Yd\X s = x s ,i,X r = x r j) can be calculated using 
the density function P{Y<i = Ud\X s = x s .i,X r = x r j) given 
in ( fTTT l and it equals h(zd). The quantities h(Yd\X s , X r ) and 
I(X s ,X r ; Y d ) have been computed and are given by (IT~8T > and 

2) Calculation of R 2 : I(X s ;Y r \X r ) is given by, 

I(X s ,Y r \Xr) = h(Yr\X r ) - h(Y r \X s ,X r ). 

The quantity h(Y r \X r be evaluated using the 

density P(Y r \X r = x r j) given in (fl4l i and h(Y r \X r ) is given 

by, 

M r -l 

h(Y r \Xr) = Yf HYr\X r = Xr,i). 

r i=0 

The quantities h(Y r \X r , X s ) = h(z r ), h(Y r \X r ) and 
I(X S ; Y r \X r ) have been computed and are given by (l20l i. (|2TT > 
and d22b . 

3) Calculation of R3: I(X s ;Y r ,Yd\X r ) is given by, 
I(X s ;Y r ,Y d \X r ) = h(Y r ,Y d \X r ) - h(Y r ,Y d \X B ,X r ), 

where h(Y r ,Yd\X r ) and h(Y r ,Yd\X s , X r ) are given by, 

M r -1 

h(Y r , Y d \X r ) = — ^2 k (Yr, Y d \X r = X r ,i)i 
r i=0 

h(Y r ,Yd\X r ,Xs) 

Af r -l*f„-l 

= WW E E h(Y r ,Y d \X r = X r ,i, Xs=Xsj). 
s r i= j= 

The quantities h(Y r , Yd\X r — x r ^,X s — x s j) and 
h(Y r ,Yd\X r = x Tj i) can be evaluated using the density func- 
tions P(Y r ,Y d \X r = x r ^,X s = x S j) and P(Y r ,Yd\X r = 
2V i) given in ( TT3T > and ( fl6] l (in the next page). 

The quantities h(Y r ,Y d \X r ), h(Y r ,Yd\X r , X s ) and 
P{X S ; Y r , Yd\X r ) have been computed and are given by d23T l. 
d24l > and d25j (in the next page). 



Appendix B 

Calculation of bounds for the HD relay channel 

WITH FINITE INPUT CONSTELLATIONS 

The probability density functions used in the calculation of 
R^-Rs are given in (I26li-(l34l> (in the next to the next page). 
Following a procedure similar to the FD relay, using (l26l - d34l >. 
i?4-i?s can be evaluated and are given by (f35T>-(f39b (in the next 
to the next page). 
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